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Abstract 

An algebraic formulation of quantum mechanics is presented. In this formulation, opera- 
tors of interest are expanded onto elements of an algebra, Q. For bound state problems in v 
dimensions the algebra Q is taken to be U{v + 1). Applications to the structure of molecules 
are presented. 


1 Introduction 

The development of new experimental techniques is giving rise to a wealth of information on com- 
plex systems. This information needs to be understood in terms of theoretical models which can, 
on one side, describe the observations and, on the other side, make predictions for other experi- 
ments. In view of the accuracy of the experiments, one needs new and more accurate mathematical 
models. For quantum mechanical problems, a natural framework is provided by the Schrodinger 
equation. In complex systems, the direct solution of the multiparticle Schrodinger equation be- 
comes rather difficult and one many seek alternative methods. In this note, an alternative method, 
called algebraic theory, will be introduced and discussed briefly. The method will then be applied 
to the study of the structure of molecules. Problems of current interest in this field are, among 
others, the substitution of atoms in a large molecule leading to a lowering of symmetry, the poly- 
merization process in which dimers, trimers, ... are formed from the original molecule, the study 
of new molecules, such as the cage fullerene molecules, C 60 , ■ • -, and, particularly important, the 
study of molecules at high excitation energy. An example of the latter will be presented. 

2 Algebraic Theory 

The logic scheme of algebraic theory is as follows: 


Quantum mechanical problem 

. U 

Algebraic Structure 


Lie Algebras 
Graded Lie Algebras 

< Infinite Dimensinal Lie Algebras (Kac-Moody) 
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f Spectra 
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In the first step, all quantum mechanical operators of interest are mapped onto the elements 
of an algebra, Q. For example, the Hamiltonian operator is written as 

H = E 0 + ^ UapGaGp + . . . , G a eQ . (1) 

Qf aft 

The algebra Q is called the spectrum generating algebra (SGA) and H is in the enveloping algebra 

otg. 

In some cases, it may happen that the Hamiltonian H contains only certain elements of Q , the 
invariant Casimir operators of Q and of one of its subalgebra chains, Q D Q f D Q n D . . . , 

H = f(C t ) . (2) 

This case, called a dynamic symmetry, plays a special role in algebraic theory, since then the 
eigenvalues of H can be obtained in closed analytic form in terms of the quantum numbers char- 
acterizing the representations of Q D G* D • • • • 

Dynamic symmetries and spectrum generating algebras have been used in various contexts for 
more than 30 years [1], As a result of the systematic investigation and use in the context of nuclear 
and molecular physics, initiated with the introduction of the interacting boson model [2] in 1974 
and of the vibron model [3] in 1981, it has become clear that all quantum mechanical problems 
in v space dimensions can be mapped onto the algebra U{y + 1) and all its states assigned to the 
totally symmetric representation [N] of U(v + 1 ) [4] . Examples of this mapping are given in the 
following section. 

3 One dimensional problems 

To begin with, consider the single case of one space dimension, v — 1. A trivial application 
of algebraic theory is provided here by the harmonic oscillator. The Schrodinger (differential) 
equation 
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( 3 ) 


with eigenvalues E n = (n + |), and eigenfunctions 

u„(x) = [ 7 ^ 2"nl] 2 ^ e~' x (4) 

(the Hermite polynomials), can be mapped into 

H=(a i a+^j , (5) 

with the same eigenvalues E m = (n -f and eigenstates 

| n >= (n!)-^(a f ) n | 0 > . (6) 



FIG. 1. The Morse potential V(r). 


Here 


The algebra is 


V2 


a = — I x + — I , a + = 


d 

dx 


(* dx) 1 


V2 

H( 2) = a, a\ 1, a + a , 


( 7 ) 

( 8 ) 
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called the Heisenberg, quantum mechanical or oscillator algebras [5]. The mapping produces a 
great simplification both in the evaluation of the matrix elements of operators which are integrals 
in the differential formulation 


Inn' — 



and algebraic functions in the algebraic formulation 

Inn' =< n I /(a, a + ) | n > 


( 9 ) 


( 10 ) 


A slightly more complicated problem is provided by the anharmonic Morse oscillator [6], Fig.l. 
The Schrodinger equation with 

H = -^-^ + V(x),V(x) = D [l-exp(-(3 x )] 2 , ( 11 ) 

with eigenfunctions 


and eigenenergies 


M*) = N.z«- v e-t+&L?- 2v - 1 {z ) 

z = 27 /e"^ , T] = ,v = 0, 1, ... ,77 - ^ , 


. „„„ [W. 1 . 1 h 2 (3 2 , i. 2 

£(«) = 2S^- ( »+-)--— (v + -) 


( 12 ) 


(13) 


can be mapped onto the algebra Q — U( 2) with elements F + , F 0 , N and corresponds to the 

dynamic symmetry U(2) D 0(2) of this algebra. The Hamiltonian is 


with eigenenergies 


H = AC , C - F 0 2 - N 2 


(14) 


E(m ) = A(m 2 — IV 2 ); m = AT, N — 2, . . . , 1 or 0 ( N = oddor even) . (15) 

The eigenvalues can be brought into the standard vibrational form introducing v = (N — m)/2 , 

N N 1 

E(y ) = — £A(Nv — u 2 ); v = 0, 1, . . . , — or — — -(N = even or odd) . (16) 

Z Z A 

With some small changes this is seen to correspond to Eq.(13). The eigenstates can be written as 
| IV, v > and observables calculated as 


< N,v' | T(G) | N,v > , (17) 

where T (G) is the appropriate operator built from the elements of Q . By making use of the 
algebraic method, all results of the anharmonic Morse oscillator can be found easily. (Note that in 
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the association of the Morse oscillator with U( 2) D 0(2 ) only the positive branch of 0(2), m > 0, 
has been used.) 

As a third simple case, consider the anharmonic Poshl-Teller oscillator 


V(x) = 


D 


cosh 2 ax 

This potential can also be associated with U( 2) j 0(2) and Hamiltonian [7] 


(18) 


H = AC 


(19) 


with eigenvalues 

E(v) = -4A(Nv - v 2 ) ( 20 ) 

One can see from (16) and (20) that the Morse and Poshl-Teller potentials have the same bound 
state spectrum (isospectral potentials). (This statement is not true in v - 2,3, ... dimensions.) 
The Morse and Poshl-Teller potentials in 1 dim belong to a class of potentials called exactly 
solvable since their eigenvalues can be written in explicit analytic form. All exactly solvable 

potentials in 1 dim have been classified. 

The overall algebraic structure of 1 dim problems can be written as 


U( 2) D 0(2) « U( 1) , 

1J. c 

H( 2) D 0(2) » 17(1) • ( 21 ) 

In this equation, the arrow with a c denotes a contraction of the algebra of U(2). (In addition 
to the contraction U( 2) c H( 2), there is another one 17(2) c E( 2), where E( 2) is the EucHdean 
algebra not discussed here.) 


4 Multidimensional problems 

In more than one dimension, the connection between the Schrodinger equation and the corre- 
sponding algebraic equation is not so straightforward, with the only exception of the harmonic 
oscillator and Coulomb problem. It is here that algebraic methods are particularly useful, since by 
formulating directly the problem in an algebraic framework one can construct the spectrum and 
calculate observables without reference to a specific form of the potential. The algebraic structure 
of three dimensional problems can be written as [4], [8] 

U( 3) D 0(3) D 0(2) (I) 

/ 


0(4) D 0(3) D 0(2) (II) 

D U( 3) D 0(3) D 0(2) (///) 


U( 4) 
^( 4 ) 
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The chain (I) corresponds to Schrodinger problems with Poschl- Teller-like potentials, the chain (II) 
corresponds to Morse-like potentials, while the chain (III) corresponds to the harmonic oscillator 
in 3 dim. 

In general, in u dim one has 

U{v) D . . . 

/ 

U{u + 1) -> ... 

\ (23) 

C (1 0(u + 1) D ... 

H(u + 1) D U{u) D ... 


where now additional chains may appear in the reduction of U(v+ 1). The five dimensional case, 
v — 5, has been extensively investigated in the context of nuclear physics [9] 


5 Algebraic Theory of Large Molecules 

The algebraic approach of Sect. 2 can be used to study molecular structure. For reasons that will 
be mentioned in the subsequent section, it is convenient to separate large molecules from small 
molecules (large here means molecules with more than 4 atoms). In large molecules each degree 
of freedom, x,y,z, is quantized with U( 2) and the total spectrum generating algebra is taken to be 
Q — Z). ©^:( 2) [10]. A calculation of spectral properties proceeds then as follows: In step 1, all 
atoms are numbered, in step 2, three coordinates are assigned to each atom for a total of 3 A/”; 
in step 3, each coordinate is quantized with U{ 2) D 0(2), thus being treated as an anharmonic 
oscillator; in step 4, the oscillators are coupled with Hamiltonians 

3jV 3j V 

H = Y h i + Y W 'J ■ (24) 

i=l i>j= 1 

The structure of the Hamiltonian (24), when written in terms of the elements of the algebra Q is 

hi — e ol + Ai(Foi — N t ) ; y>ij = \ij(F + iF_j + F_ , F + } ) (25) 

The hi terms are diagonal in the basis t/,( 2) D 0,(2) characterized by the quantum numbers 
I Ni, v i > discussed in Sect. 3, while the w.j are given in Ref. [10], [11]. In the final step 5, the 
spurious species corresponding to overall rotations and vibrations, are identified and removed by 
diagonalizing the Hamiltonian 


H' = H + XV , (26) 

where F is a projection operator into the spurious species and A is taken to be a large number 
such that the spurious species are moved to a large energy. The removal of the spurious species 
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leaves 3Af — 6 non spurious vibrations. This procedure produces the vibrational spectrum of the 
molecule. 


H 



FIG. 2. Schematic representation of benzene ( C§H§ ). 

In a similar way one can compute intensities of transitions. There are two types of transitions 
of importance in molecules, infrared (IR) and Raman (R) transitions. For infrared transitions, 
the appropriate operator is a vector. Each component x,y,z of this vector ts written in terms of 
elements of Q, i.e. 


3A r 

T x = ^2 a itX U 

i=i 


l _ e -0i(F +i +F_i) 


(27) 


The matrix elements of the operator (27) (and T y , T z ) are then evaluated algebraically. For 
Raman transitions, the appropriate operator is a symmetric quadrupole tensor. The six compo- 
nents of this tensor, a;®, y 2 , z 2 , xy, xz, yz , are also written in terms of Q and their matrix elements 
evaluated algebraically. 

As an example of vibrational analysis of large molecules, consider the case of benzene, C 6 H 6 . 
The benzene molecule has the geometric structure shown in Fig. 2. A problem which arises in 
large molecules is that of the discrete symmetry of the molecule. In the case of benzene, the 
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appropriate symmetry is V 6 /,. The discrete symmetries of molecules can be simply implemented 
in the algebraic framework. For example, consider the six stretching vibrations of the hydrogen 
atoms in benzene. All hydrogen atoms are equivalent. The Hamiltonian (24) which describes 
those vibrations g 

H = Y t A i Ci+ £ , (28) 

»=i i<j= i 

where C, and Mij are a short-hand notation for the terms in (25), must be such that one cannot 
distinguish the equivalent atoms. Thus, all A,’s must be equal, A t = A. In the interaction term, 
there are three contributions, first, second and third neighbor interactions. These too must be 
equal, A^ = A^, A t -j^ = A and A-j 7 ^ = \( in h The Hamiltonian H thus becomes 

H = AC+X^S^ + A + A( J// ^( //J ) , (29) 

and is characterized by a smaller set of parameters. The operators are ap- 

propriate linear combinations of the CV s and Af^’s. A corollary of the algebraic method is that 
certain linear combinations of the operators Af,j are symmetry adapter operators of T> 6 h and the 
irreducible representations of V 6 h are obtained automatically by diagonalizing them [10], [12]. 

Using the algebraic method discussed above it has been possible to study the complete spec- 
troscopy of benzene [13]. This molecule has 12 atoms and thus 36 degrees of freedom, 6 of which 
are spurious. The 30 non spurious species are shown in Table I. 


TABLE I: Coordinates and symmetry species of benzene. 


Coordinates 

Number 

Species 

Spurious 

CH strekch 
CC stretch 
CH in plane bend 
CH out of plane bend 
CC in plane bend 
CC out of plane bend 

6 

6 

6 

6 

3+3 

3+3 

E 2g -f B\ u + E\ u -f A\ g 
A\g + B 2u + El u + E'2g 
E\u + B 2 U ~t~ E 2 g + A 2 g 
A 2u + B 2g + E\g + E 2u 
B\ U + E 2 g 

B 2g + E 2g 

i 

A\g + E\ u 
A\ n + E\g 


t 


The calculation describes the observed vibrational states not only in the low excitation energy 
region, fundamental vibrations, but also in the high excitation energy region, overtones. An 
example is shown in Fig. 3. This region cannot be described in the harmonic approximation and 
thus the use of algebraic methods based on the anharmonic Morse or Poschl- Teller oscillators are 
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crucial for an accurate description of the observed spectra. 



8725 8775 8825 8875 

WAVENUMBERS [cm -1 ] 


FIG. 3. (a) Opto- thermal spectrum in the region of the Av = 3 overtone of the 
stretching CH mode of benzene. The full-line is a low-resolution experiment, (b) The 
spectrum calculated by means of algebraic theory. From R.ef.[14]. 
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6 Algebraic Theory of Small Molecules 

For molecules with a number of atoms less or equal to four, it is possible to quantize each vector 
degree of freedom, r = (x,y,z), directly in terms of 17(4). When quantized in this way both 
rotations and vibrations are simultaneously included. 



FIG. 4. Bond variables for small molecules. 


It is also convenient to treat as vector variables the bond degrees of freedom, Fig.4, rather than the 
coordinates of the individual atoms, thus avoiding the problem of spurious states. The quantization 
scheme is thus here Q = X), ©17(4). This scheme has been extensively used to treat diatomic 
molecules with 17(4), triatomic molecules with 17(4) © 17(4) [15] and four atomic molecules with 
17(4) © 17(4) © 17(4) [16]. It has also been possible recently to study high order interactions such 
as rotation- vibration couplings. 

7 Conclusions 

Algebraic methods have been used in recent years in the study of molecular structure. When 
applied to this system, algebraic theory offers two main advantages: (i) The use of Lie algebras to 
describe the interaction (Morse, Pos chi- Teller, . . .) allows one to extend the traditional harmonic 
analysis to anharmonic analysis. One can thus deal easily with highly excited states of molecules 
where anharmonicities play a crucial role (a subject of current experimental interest especially in 
connection with intramolecular relaxation and energy transfer.) (ii) The use of algebraic operators 
to couple the individual modes of a molecule allows one to construct symmetry adapted states in 
a simple way. One can thus deal with complex molecules where discrete symmetries play a crucial 
role. 

The algebraic method can be used in molecules in two ways: 

A) With rotation and vibrations treated separately. In this case the spectrum generating algebra 
is 
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G = Grot © Gvib 

GviB = £ ®I/t(2) • 


( 30 ) 


1 


B) With rotations and vibrations treated simultaneously. In this case the spectrum generating 
algebra is 


£ = £©£/, (4) • (31) 

The latter case is more complete, but more difficult to treat than the former, since one has to deal 
with the Racah algebra of U(4). 

In view of its simplicity, the method is particularly well suited for a study of complex systems 
such as macromolecuhjs, clusters, polymers, . . .. Work in this direction is in progress. An account 
of the algebraic theory of molecules is given in Ref.[17] and the mathematical formalism of Sects. 
5 and 6 is reviewed in Ref. [18]. 
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